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This report provides an overview of the Hyper-X research vehicle Monte Carlo analysis 
conducted with the six-degree-of-freedom simulation. The methodology and model 
uncertainties used for the Monte Carlo analysis are presented as permitted. In addition, the 
process used to select hardware validation test cases from the Monte Carlo data is described. 
The preflight Monte Carlo analysis indicated that the X-43A control system was robust to 
the preflight uncertainties and provided the Hyper-X project an important indication that 
the vehicle would likely be successful in accomplishing the mission objectives. The X-43A in- 
flight performance is compared to the preflight Monte Carlo predictions and shown to 
exceed the Monte Carlo bounds in several instances. Possible modeling shortfalls are 
presented that may account for these discrepancies. The flight control laws and guidance 
algorithms were robust enough as a result of the preflight Monte Carlo analysis that the 
unexpected in-flight performance did not have undue consequences. Modeling and Monte 
Carlo analysis lessons learned are presented. 


Nomenclature 


alpha 

= angle of attack 

AIL 

= aircraft- in-the- loop 

AMW 

= all-moving wing 

Avg 

= average 

CA 

= axial force coefficient 

CA bias 

= axial force bias 

CFD 

= computational fluid dynamics 

CG 

= center of gravity 

Cl 

= roll moment coefficient 

Cl bias 

= roll moment bias 

Clp 

= roll moment caused by roll rate 

Cm 

= pitch moment coefficient 

Cm bias 

= pitch moment bias 

Cmo 

= pitching moment at zero degree angle of attack 
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Cmq 

= 

pitch moment caused by pitch rate 

Cn 

= 

yaw moment coefficient 

Cn bias 

= 

yaw moment coefficient bias 

Cnp 

= 

yaw moment caused by roll rate 

Cnr 

= 

yaw moment caused by yaw rate 

CP 

= 

pressure coefficient 

CY 

= 

side force coefficient 

CY bias 

= 

side force bias 

CZ 

= 

normal force coefficient 

CZ bias 

= 

normal force bias 

DFRC 

= 

Dryden Flight Research Center 

FADS 

= 

flush airdata sensing 

FMU 

= 

flight management unit 

G 

= 

gravitational acceleration 

GNC 

= 

guidance, navigation, and control 

GRAM 

= 

Global Reference Atmospheric Model 

h 

= 

altitude 

HIL 

= 

hardware-in-the-loop 

HXLV 

= 

Hyper-X Launch Vehicle 

HXRV 

= 

Hyper-X Research Vehicle 

INS 

= 

inertial navigation system 

kft 

= 

thousands of feet 

MCAT 

= 

Monte Carlo analysis tool 

PID 

= 

parameter identification 

PPT 

= 

pressure port transducer 

SepSim 

= 

detailed simulation of separation dynamics 

SIM 

= 

simulation 

UA 

= 

unavailable 

6-DOF 

= 

six degree-of-freedom 

G 

= 

standard deviation 


I. Introduction 

B efore the initial flight-testing of a new aircraft, there is always some uncertainty as to how the vehicle will 
actually perform. These uncertainties increase when the vehicle is a unique aircraft flying in a seldom-explored 
flight regime. One such unique vehicle is the X-43A, which was the first unmanned vehicle to fly at hypersonic 
speeds with a scramjet engine. A Monte Carlo analysis was conducted to predict and bound the performance of the 
X-43A, hereafter referred to as the Hyper-X research vehicle (HXRV). Monte Carlo analysis is a useful and 
commonly used technique to assess system performance to predicted parameter variations by means of a stochastic 
process. The analysis consists of independent and random variations of all applicable model inputs and parameters 
within the estimated uncertainty bounds. Random combinations of input uncertainties can produce different, and 
potentially worse, results than deterministically varying a single uncertainty or set of uncertainties. Monte Carlo 
analysis consists of making multiple simulation runs, with each run containing a unique set of parameter variations, 
and analyzing the results. A close examination of these simulation runs provides insight into the system performance 
and characteristics. Monte Carlo analysis forces the flight control laws and guidance algorithms to be robust to a 
large range of uncertainties. Stress testing past the modeled uncertainty range can provide an assessment of the 
vehicle robustness to the unknown unknowns present in every flight research program. This type of analysis is 
particularly useful for unmanned vehicles where there is no pilot in the loop acting as an adaptive controller. 

The purpose of this report is to provide an overview of the HXRV Monte Carlo analysis tools, present the 
uncertainties used as part of the Monte Carlo analysis (as permitted), and discuss the results of the Monte Carlo 
analysis. In addition, limited comparisons between simulation and flight data are presented with an emphasis placed 
on highlighting modeling shortfalls. The selection process for the subset of runs used during hardware validation 
testing is described, and lessons learned and observations concerning the Monte Carlo analysis are presented. 
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II. Project Description 

The goal of the NASA Hyper-X program was to demonstrate an advanced, airframe-integrated, airbreathing, 
hypersonic propulsion system in flight, thus validating the supporting tools and technologies. 1 The program sought 
to design, build, and flight-test a series of three small, autonomous, scramjet-powered Hyper-X research vehicles 
(HXRV), designated X-43A, at Mach 7 and 10. The vehicles were dropped from the NASA Dryden Flight Research 
Center (DFRC) B-52B airplane (The Boeing Company, Chicago, Illinois) over the Pacific Ocean and rocket-boosted 
to the test point. Prior to flight test, the Hyper-X program conducted extensive ground testing, simulation, and 
analysis. In support of these activities, a high-fidelity six degree-of-freedom (6-DOF) nonlinear simulation of the 
research vehicle was developed at DFRC. The research vehicle simulation was meant to accurately model the 
vehicle dynamics after separation from the Hyper-X launch vehicle (HXLV), including the engine test and 
controlled flight down to the water. Many of the simulation models included uncertainties or dispersions that were 
used in the Monte Carlo analysis. 

A. The Hyper-X Research Vehicle (HXRV) Description 

As shown in Fig. 1, the HXRV was an unmanned autonomous vehicle that measured approximately 12 ft long 
and 5 ft wide, and weighed approximately 3000 lb. The scramjet engine was attached to the underside of the HXRV 
and contained a cowl door on the leading edge of the engine that was used to control airflow through the engine. The 
HXRV had four control surface effectors: a left and right all-moving wing (AMW) and twin rudders. The HXRV 
was fitted to the front end of the Hyper-X launch vehicle (HXLV), a modified Pegasus 8 (Orbital Sciences 
Corporation, Chandler, Arizona) rocket used to boost the HXRV to the desired test condition. The combination of 
the HXRV and HXLV was known as the X-43A stack, and was carried to the launch point under the wing of the 
NASA DFRC NB-52B airplane, ship number 008. 

1. Flight Test of the Hyper-X Research Vehicle 

Three HXRVs were built for the Hyper-X program. Two of the vehicles were intended for missions at Mach 7, 
and one vehicle was intended for a mission at Mach 10. The HXRVs collectively received the designation of X-43A. 
All three vehicles have the same outer mold line; the primary difference among the vehicles is the internal engine 
flowpaths. Flight 1 was intended to reach Mach 7 and was attempted on June 2, 2001. The HXLV lost control 
shortly after launch, resulting in the loss of both the HXLV and HXRV. 2 During Flight 2, flown to Mach 6.8 on 
March 27, 2004, the HXRV successfully demonstrated the in-flight operation of the scramjet. 3 All of the goals for 
that mission were achieved, including positive acceleration of the vehicle by the scramjet. The third and final flight 
was flown to Mach 9.6 on November 16, 2004. 4 * During both successful missions, the HXRV was in a controlled 
autonomous flight from the point of separation to splashdown in the Pacific Ocean. 

2. The Hyper-X Research Vehicle Mission Description 

The HXRV mission timelines were similar for Flight 2 and Flight 3. The principal difference between the two 
flights is the higher Mach number attained during Flight 3, which resulted in a correspondingly longer descent 
trajectory. Figure 2 shows an overview of the HXRV mission. The launch, scramjet experiment, and descent 
portions of the mission were conducted off the coast of southern California. For both missions, the HXLV propelled 
the X-43A stack to a separation altitude of approximately 100,000 ft and a dynamic pressure of approximately 
1000 psf. The HXRV separated at Mach 6.9 during Flight 2 and at Mach 9.7 during Flight 3. Several seconds after 
separation, the primary experiment, which consisted of the scramjet engine test, was conducted. The primary 
experiment lasted approximately 45 s for Flight 2 and 30 s for Flight 3. After the completion of the primary 
experiment, the cowl door was closed and a recovery maneuver was initiated to arrest the vehicle descent and 
dynamic pressure buildup. After the descent was arrested, an unpowered trajectory was flown to a splashdown into 
the Pacific Ocean. During the descent, a series of parameter identification (PID) maneuvers were performed at every 
integer Mach number down to Mach 2 , 5,6 These maneuvers began at Mach 5 for Flight 2 and at Mach 8 for Flight 3. 

III. Monte Carlo Analysis Overview 

The DFRC 6-DOF nonlinear simulation of the HXRV was used in conjunction with a series of Matlab 8 (The 
MathWorks, Natick, Massachusetts) scripts to conduct the Monte Carlo analysis. The primary objectives of the 
Monte Carlo uncertainty analysis were to assess vehicle stability in the presence of modeling and measurement 
uncertainties, the flight conditions during the engine test, the unpowered descent trajectory, and the controlled flight 
impact area. Statistical information from the Monte Carlo analysis was examined, but was not the primary product. 
The project was most concerned with assessing the vehicle performance and the expected flight envelope in the 
presence of the expected model uncertainties. This performance assessment turned out to be largely Boolean in 
nature. The vehicle performance was either acceptable or not. Software performance, including the flight control 
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laws, guidance, and scramjet controller 7 was a part of this assessment. The Monte Carlo objectives were achieved by 
conducting many thousands of simulation runs, with each run having a random combination of modeling 
uncertainties. To gain a greater understanding of the relationships between the dispersions and their effects on the 
trajectory, data were correlated to determine combinations of factors that most influenced the performance of the 
aircraft. A secondary objective consisted of using the uncertainty analysis data to find a limited set of runs that could 
be used for thorough validation testing of the flight hardware and software systems. 

A. The Hyper-X Research Vehicle Monte Carlo Analysis Tool (MCAT) 

A Monte Carlo analysis tool (MCAT) was developed to stress the vehicle system with predicted uncertainties. 
The HXRV MCAT tests the vehicle performance in the presence of random and unique sets of input parameter 
uncertainties. The MCAT consists of a set of Matlab 17 M-files that are used as controlling scripts to drive the 
nonlinear simulation with various dispersions, analyze the results, save the data, and plot the data for further 
analysis. Matlab® M-files are used to automate the process of generating the thousands of simulation scripts. 

An execution of the MCAT consists of a large number of simulation runs with the first run being nominal and 
the rest containing different combinations of uncertainties. 

Manipulation and storage of the Monte Carlo results presented numerous challenges. Several gigabytes of data 
were produced by each set of simulation runs. The analysis tools run in the Matlab® workspace were incapable of 
handling a complete Monte Carlo dataset. The large amount of data resulted in the partitioning of the Monte Carlo 
data into multiple .mat files. A separate .mat file was generated for each simulation output variable and this file 
contained data from all of the simulation runs. To minimize the amount of data stored from a Monte Carlo analysis, 
trade-offs were made between the rate at which the data were saved and the number of signals saved for short-term 
and long-term storage. The parameters of greatest interest were saved at 100 Hz during the engine experiment phase. 
During other phases of the simulated missions, data were saved at 10 Hz. Relatively slowly-varying parameters, 
such as altitude, Mach, latitude, and longitude were stored at 1 Hz. Table 1 shows the Monte Carlo parameters saved 
at each rate. Parameters of greatest interest were retained indefinitely, whereas less important data was deleted after 
a week to allow storage space for the next set of Monte Carlo results. 

Figure 3 is a flowchart graphically representing the stochastic MCAT process. The inputs of each simulation run, 
including the desired uncertainties, are defined in an input file that is called a simulation input script. The MCAT 
generates the simulation input scripts and details the number of simulation runs. The MCAT then uses the 
simulation input scripts to run the simulation in a batch mode a number of times specified by the MCAT and saves a 
subset of the simulation parameters that are recorded for later analysis. The simulation input scripts were deleted 
after each Monte Carlo run to save space. The simulation output data of each of these saved parameters are collected 
from each simulation run into separate matrices (one for each parameter). The data matrices are then stored for 
future analysis and plotting. 

For the MCAT, every input parameter that is varied has an uncertainty expressed in terms of a probability 
density function. The probability density functions for each input parameter are expressed in terms of its distribution 
type (uniform or Gaussian), a mean, and a standard deviation. The vast majority of the simulation input parameters 
are varied independently. Several parameters, such as the moments of inertia, were checked for consistency. If a 
parameter was found to be inconsistent with the other parameters with which it is correlated, the value was thrown 
out and recast until a consistent value was found. The input parameter values are based upon their particular 
probability density functions. They are selected within the ± 3 sigma bounds for each run with either a uniform or 
Gaussian distribution (the distributions are specified by the user in the MCAT). If a distribution is Gaussian, the 
3 sigma bounds encompass 99.87 percent of the possible values of the uncertain parameter. The uncertainties 
associated with a parameter were provided by each discipline responsible for that portion of the model (i.e. 
aerodynamics, propulsion, flight control, etc.). 7 

As a part of each Monte Carlo analysis, summary statistics were extracted from the dataset and placed in a 
Monte Carlo report. A temporary 100 Hz dataset was created to extract the statistical information from the highest 
simulation data rate available. The temporary dataset was deleted after the statistics had been gathered. A number of 
statistics, including the maximum and minimum values, were recorded for certain parameters during particular flight 
phases. In addition, the run number associated with each maxima or minima was recorded. Table 2 shows a sample 
page from the Monte Carlo report. This report was examined and compared to previous reports to determine if any 
unusual cases or events had been generated in the latest set of Monte Carlo runs. Effective use of these summary 
statistics helped to reduce the long-term data storage requirements and identify runs of interest. The simulation 
scripts for each Monte Carlo analysis could be easily regenerated. This allowed for a simulation run to be repeated 
and more data acquired if, at a later date, a review of the summary or lower rate data indicated an item of interest. 
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In preparation for the second Mach 7 flight attempt, 1600 Monte Carlo simulation runs were conducted with the 
latest models and software as part of the final analysis. Many thousands of Monte Carlo simulation runs had been 
conducted prior to the generation of the final dataset as the analysis had been run continuously as the simulation 
models had been refined and the vehicle software matured. The number of simulation runs used for the final analysis 
was determined by examining the statistics of notable parameters, such as the maximum value for filtered normal 
acceleration, seen in Fig. 4. The asymptotic trend of the parameters showed that there would be little expected 
change in the output of the Monte Carlo analysis when performing more than 1600 runs. 

IV. Simulation Description 

The Hyper-X program developed a high-fidelity 6-DOF nonlinear simulation model of the HXRV. The HXRV 
simulation does not model the HXLV and was usually initialized to the separation conditions from the HXLV. A 
separate 6-DOF simulation was developed to model the HXLV, and is beyond the scope of this report. A simplified 
model of the separation dynamics is included in the HXRV simulation. A separate, more detailed simulation, known 
as SepSim, was developed to fully model the HXRV separation from the HXLV and is beyond the scope of this 
report. The HXRV simulation was based on the standard DFRC simulation architecture 8 and contained detailed 
models that included the vehicle aerodynamics, mass properties, atmospheric, and wind models. Vehicle subsystems 
modeled in the simulation included the flight control system, actuator, engine and fuel systems, the inertial 
navigation system (INS), and the flush airdata sensing (FADS) system. The flight management unit (FMU) software 
models consisted of guidance, navigation, and flight control systems, in addition to the propulsion system control. 
Since most of the components in the research vehicle were single-string (not fault redundant), no failure modes were 
modeled in the simulation. The majority of the simulation core source code and models are written in Fortran 77 
with several of the models consisting of autocode generated from Simulink ' (The MathWorks, Natick, 
Massachusetts) models. The simulation graphical user interface is written in C programming language. 

The simulation can be run in either a real-time mode or a non-real-time batch mode. The batch mode operation is 
faster than real-time and was used for most engineering analysis including Monte Carlo analysis. The real-time 
mode was used only when the simulation was connected to ground test hardware. The base simulation frame rate is 
200 Hz. Vehicle models and simulation routines are called at the appropriate multiples or submultiples of the base 
frame rate. Numerical integration of the simulation state information is performed using a second-order fixed-step 
Runge-Kutta algorithm optimized for fast execution. Oblate-Earth, 6-DOF equations of motion are employed. The 
simulation employs a command line interface along with a scripting capability. The scripting capability allowed for 
the simulation execution of script files containing multiple commands that had been strung together. The MCAT 
made extensive use of the simulation scripting capability to generate a unique simulation setup script for each Monte 
Carlo simulation run. Nearly every variable in the simulation is capable of being recorded in an external data file. 
The selection of output signals was done through the script files. 

V. Model Uncertainties 

The Monte Carlo analysis consisted of applying uncertainties to the simulation vehicle models. A total of 
286 model uncertainties were used as a part of the final analysis for the second and third flights of the HXRV. The 
uncertainty distributions for each model were estimated using measurement error buildup, test data spread, or 
engineering judgment. The model uncertainties were generally developed by the engineering discipline team 
responsible for producing the model. These model uncertainties could be grouped in the following categories: 
separation condition, mass properties, control surface actuator performance, HXRV aerodynamics, separation 
aerodynamics, separation piston performance, FADS, atmosphere, winds, HXRV sensors, engine performance, and 
propulsion system performance. Uniform uncertainty distributions, described by their minimum and maximum 
values, were generally used as they produced the most conservative results. Uniform distributions were used for 
separation conditions, engine performance, winds, and inertial navigation system (INS) errors. To describe each 
Gaussian distribution, four numbers were used: mean, standard deviation (1-sigma), minimum, and maximum. The 
minimum and maximum numbers were used to limit the dispersions to better match their intended distributions and 
to maintain a 3 -sigma limit on the applied uncertainties. The Matlab 1 ' rand or randn function was used to generate a 
random number with either a uniform or Gaussian distribution for each uncertainty. These random numbers were 
then linearly mapped to the model uncertainty distribution to generate a model uncertainty parameter for input to the 
simulation. 

An overview of the Monte Carlo uncertainties is presented in Table 3. A detailed discussion of each area 
follows, with details presented where the data is not proprietary, sensitive, or otherwise restricted in distribution. 
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The HXRV simulation was initialized at the point of separation from the HXLV for the Monte Carlo analysis. 
The flight conditions used to initialize the HXRV simulation were derived from the project requirements document 
and the HXLV Monte Carlo simulation results. The HXLV Monte Carlo results generally fell within specified 
bounds. When a specific project requirement existed, that bound was used to determine the mean and ± 3 sigma 
bounds for a particular condition since the requirement was more conservative than the HXLV Monte Carlo results. 
When no preexisting project requirement existed, the output of the HXLV Monte Carlo simulation was used. 
Table 4 details the separation conditions and uncertainties used to initialize the HXRV simulation for the Monte 
Carlo analysis. 

The mass properties uncertainty model was developed by the HXRV structures team and is largely based on the 
accuracy of the test data. The vehicle weight, center of gravity, and principle moments of inertia were measured on 
the ground and adjusted for the final flight configuration. Measurement uncertainties and engineering judgment were 
used to determine the appropriate uncertainties. The mass property uncertainty values are presented in Table 4. 

The HXRV control surface actuator models are discussed in detail in Ref. 9. A high- fidelity, proprietary model 
of the HXRV control surface actuators was provided by the manufacturer, Moog (New York) and updated by the 
HXRV guidance, navigation, and control (GNC) team, which also developed an uncertainty model for the actuator 
model. This high-fidelity model was used for the HXRV Monte Carlo analysis for Flights 2 and 3. The HXRV GNC 
team also developed a simplified second-order model that approximated the actuator performance. This simplified 
model was used for Monte Carlo analysis prior to the first flight attempt and these model uncertainties are presented 
in Table 4. In general, the simplified actuator uncertainty model is more conservative than the high-fidelity model. 

The HXRV vehicle aerodynamics model was developed from wind tunnel testing and computational fluid 
dynamics (CFD). 10 The aerodynamic uncertainty model was developed based on the wind tunnel testing and CFD 
results." In addition, work done in developing the X-33 uncertainty model 12 was leveraged in the uncertainty 
development. The aerodynamic coefficient uncertainty is a function of Mach number and engine operation. The 
aerodynamic uncertainties for cowl-closed flight decrease as the Mach number decreases. The aerodynamic 
uncertainties are increased for cowl-open flight and increased further while the scramjet is operating. The 
uncertainties in the principle aerodynamic coefficients during cowl-closed flight are presented as percents in 
Table 4. Aerodynamic coefficient biases were used, but are not presented as they were not formulated in terms of a 
percent. The magnitude of the aerodynamic coefficient biases is generally small relative to the coefficients. 

A simplified separation aerodynamics model was developed which applied delta force and moment increments 
to the free-flight HXRV aerodynamic coefficients to approximate the aerodynamic interactions expected to occur 
during the separation event. This model is known as the Sep Delta model. Uncertainties on these delta increments 
were developed and incorporated into the Monte Carlo analysis. The Sep Delta model was derived from SepSim 
analysis. A detailed presentation of the Sep Delta model along with its uncertainties is beyond the scope of this 
report. The principle purpose of the Sep Delta model in the HXRV Monte Carlo simulation was to replicate the 
spread in flight conditions the vehicle could be expected to encounter after the separation event. A detailed 
discussion of the SepSim analysis and results is contained in Ref. 13. 

Two pistons were used to push the HXRV away from the HXLV during the separation event. A piston force 
model was developed, along with uncertainties, primarily from test data. A piston performance model was 
implemented in the HXRV simulation, in addition to SepSim. A detailed discussion of the piston performance 
model and uncertainties is beyond the scope of this report. 

A flush airdata sensing (FADS) system was employed on the HXRV. 14 The FADS system provided an angle of 
attack estimate that was blended with the inertial measurement following the engine test during the Mach 7 flight. 
The FADS system was only used real-time during the unpowered descent portion of the Mach 7 flight. The FADS 
performance uncertainties were examined by introducing uncertainties in the simulated pressure port measurements. 
There are three components to the FADS pressure port uncertainties. White noise was added to each frame as a 
uniformly-distributed number continuously varying from -1 to 1. There was also an uncertainty in the pressure 
coefficient at each port, which was a function of Mach and alpha. Finally, a PPT (pressure port transducer) 
measurement uncertainty was added as a normally-distributed number between -1 and 1. For the Monte Carlo 
analysis, uncertainties were applied only to the pressure ports used as part of the angle of attack estimation 
algorithm. Table 4 details the FADS uncertainties. 

The HXRV simulation contains a separate atmosphere model for each month. This model was developed by the 
DFRC meteorological group and tailored for the expected HXRV engine test location. The model is derived from 
Global Reference Atmospheric Model (GRAM-95) data. 15 The atmosphere model used by the HXRV onboard 
software is constant during a Monte Carlo run and is set to the expected flight month. Atmospheric uncertainties are 
modeled by choosing a uniform random month for the atmosphere model in the nonlinear simulation. By randomly 
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varying the simulation atmosphere model between the 12 available models, adequate uncertainty was applied to the 
HXRV onboard software relative to the simulated values. 

North-South, East-West, and vertical wind components are independently varied every 10,000 ft. North-South 
and East-West wind data was provided by the DFRC meteorology group. Vertical wind mean data is from 
GRAM-95, and also varies by altitude and month. Because of sparse statistical data for vertical winds over the 
ocean, a worst-case 3 sigma value of 10 ft/s was used. Table 4 provides sample wind uncertainties. 

A high-fidelity sensor model, including uncertainties, was developed to model the performance of the INS 
contained in the H-764 Honeywell flight control computer (Honeywell, Inc., Minneapolis, Minnesota). The INS 
package is capable of utilizing global positioning system (GPS) updates to aid the position and attitude solution, but 
the HXRV used the pure inertial solution throughout its mission, so GPS performance was not included in the sensor 
model. Information used to develop the high-fidelity sensor model is proprietary and therefore, the model details are 
not presented here. 

The HXRV simulation contains an engine database that approximates the steady-state scramjet engine thrust and 
moments, hereafter referred to as performance. 16 The simulation engine database consists of maximum and 
minimum values as a function of flight condition and fuel-to-air ratio. Engine performance is varied for each run by 
setting an engine performance parameter that selects a point between the minimum and maximum values in the 
database. Uniformly distributed values between 0 and 1 were used to vary the engine performance. The engine 
performance parameter was set at the beginning of a Monte Carlo simulation run and not varied thereafter. A value 
of 0.5 corresponds to the midpoint of the maximum and minimum values and is considered the nominal engine 
performance. 

Uncertainties in the HXRV scramjet fuel system performance were modeled in the Monte Carlo analysis. 
Uncertainty in the fuel flow, igniter flow, and airflow were also modeled. 

Most of the applied uncertainties were independent of each other; however, some reasonableness checks were 
performed on correlated parameters. These included angle of attack, pitch attitude, and flight path angle, as well as 
altitude, velocity, and dynamic pressure. Moments of inertia were also evaluated to ensure that the sum of two 
principle axes would be greater than the third axis. 

In addition, the HXRV simulation random number seed was set by the Monte Carlo tool. The random number 
seed is used to initialize the random functions in the HXRV simulation that apply noise to the sensors and FADS 
pressure measurements. This explicit setting of the random number seed allowed for repeatability of a single Monte 
Carlo run. 


VI. Mission and Trajectory Sensitivity Analysis 

Key relationships between uncertainties and performance were determined using a sensitivity study. In general, 
the project was most interested in determining the uncertainty combinations that would produce the maximum or 
minimum of a particular metric. For example, one metric was the vehicle positive acceleration caused by the 
scramjet engine operation. As a result, the maximum and minimum positive acceleration values were of interest 
along with any key uncertainties influencing the value. 

Performance metrics were unique to the various mission segments. The following were some of the metrics that 
were tracked during the separation event and parameter identification maneuvers: surface command rates and 
positions, vehicle attitude rates and angles, and normal acceleration. The most important performance metrics were 
the enginc-on flight conditions: primarily angle of attack, Mach number, and dynamic pressure. These metrics were 
used to provide the scramjet team with the expected flight conditions for their own detailed, off-nominal analysis of 
the scramjet performance. The robustness of the guidance algorithm to off-nominal engagements was assessed by 
examining the impact of off-nominal headings when the algorithm was engaged. Over the rest of the descent to the 
ocean, surface commands and positions, altitude error, vehicle attitude, and flight condition were tracked. The 
vehicle orientation, velocity, splashdown point, aim point error, and time aloft were then determined at the end of 
each run. This final information was used to evaluate the performance of the guidance routine and to develop a 
splash point boundary for range clearance. 

Maximums and minimums of these parameters of interest were drawn on scatter plots versus each model 
uncertainty input for a set of Monte Carlo runs to determine the uncertainty parameters to which the vehicle 
performance was most sensitive. An obvious example is seen in Fig. 5, which shows how increasing the drag (axial 
force aerodynamic coefficient) decreases the distance traveled. 

The HXRV was found to be sensitive to a model uncertainty if it caused the vehicle to reach a maximum or 
minimum limit of some parameter of interest, such as angle of attack, normal acceleration, or control surface 
deflection. In addition, if the vehicle exhibited an undesirable behavior, such as a pitch or roll oscillation, the 
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cause(s) were noted as parameters of interest. Before the first flight attempt, a ranked list was formed of the model 
uncertainties to which the HXRV performance metrics and trajectory were most sensitive. 

1. Separation angle-of-attack 

2. Separation sideslip angle 

3. Separation pitch attitude 

4. Separation heading 

5. Separation altitude 

6. Separation Mach 

7. Lateral center of gravity location 

8. Actuator damping ratio 

9. Axial force aerodynamic coefficient (approximately equivalent to Drag) 

10. Normal force aerodynamic coefficient (approximately equivalent to Lift) 

11. Pitching moment aerodynamic coefficient 

12. Error in pitch attitude caused by mounting and measurement 

Stress testing was conducted with combinations of these parameters being set at or beyond their minimum or 
maximum expected values. This testing resulted in a handful of cases that stressed the vehicle performance at 
varying points throughout the mission. These limited numbers of cases could be run quicker than a full Monte Carlo 
analysis to assess vehicle performance, while still stressing the vehicle software and performance. These stress tests 
were used to quickly examine design changes to the flight control laws and guidance algorithm. In addition, the 
stress test analysis partially accounted for the unknown unknowns present in every flight research program by 
examining the vehicle performance with what are believed to be unrealistic model uncertainties. These tests 
provided an assessment of the vehicle robustness in excess of that which is expected to be required. 

VII. Validation Testing 

Two types of flight hardware and software validation tests were performed: hardware-in-the-loop (HIL) and 
aircraft- in-the- loop (AIL) tests. 17 The purpose of HIL testing was to validate the hardware and software performance 
of the flight control computer, which is a portion of the flight management unit (FMU). For HIL testing, the FMU 
was connected to the simulation. Simulated inertial sensor data was provided to the FMU by the simulation. The 
flight software running on the FMU generated control surface commands and propulsion control valve commands 
that were in turn fed to the simulation. The number of test cases run in the HIL environment was limited because of 
the test complexity and time involved. The AIL tests were similar to the HIL tests but also included the vehicle 
flight actuator controller and actuators as part of the test configuration. The AIL testing placed flight hardware at 
risk, so an even more limited number of test cases were chosen for this type of testing. 

Taking these factors into consideration, a small set of runs was chosen for validation testing from the thousands 
of Monte Carlo runs. Approximately 20 runs were performed using the HIL setup, and less than 5 runs were 
conducted as a part of AIL testing. The runs chosen for HIL and AIL testing were chosen from a Monte Carlo 
dataset by finding the minimum number of runs that bounded some desired aspect of the mission. For example, 
Fig. 6 shows the 5 cases found that bounded 90 percent of the altitude profiles from 1600 Monte Carlo cases. In 
general, the HIL cases were chosen to bound the flight conditions seen in the Monte Carlo analysis and the AIL 
cases were chosen to stress some aspect of the control surface actuators such as response to a large step change, 
maximum/minimum deflection, or high levels of surface motion. An automated tool was developed that selected the 
minimum number of cases that met the desired criteria. Table 5 shows the criteria the tool maximized or minimized 
during a given mission segment. The tool worked its way through the Monte Carlo dataset and selected a run for 
HIL testing if it was 90 percent of the way to a specified bound. Cases exceeding the 90 percent threshold for 
multiple specifications were ranked higher and the highest ranked cases were chosen for HIL analysis. 

VIII. Flight Data Comparisons to Preflight Monte Carlo Results 

A summary of the HXRV in-flight performance comparison with preflight Monte Carlo results is discussed for 
the Mach 7 and Mach 10 flights, with an emphasis placed on highlighting modeling shortfalls. Reference 18 
provides a detailed discussion of the HXRV performance and comparisons with Monte Carlo data during the Mach 7 
flight and Ref. 19 provides the same for the Mach 10 flight. 


8 


American Institute of Aeronautics and Astronautics 

092407 



Flight results are presented for the time period after the separation event. The HXRV simulation only roughly 
approximates the separation event and this time period is fully modeled in detail by the HXRV separation 
simulation. 

A. The Mach 7 Mission 

For the Mach 7 flight, the in-flight performance, when compared to Monte Carlo predictions, indicates that the 
analysis was fairly accurate at predicting actual flight performance and bounds. Notable differences between the 
Monte Carlo predictions and the in-flight performance are discussed. Figure 7 displays the inertial angle of attack 
during the scramjet test. In Figs. 7-15, flight data is compared to a nominal simulation run started at the in-flight 
separation point. Monte Carlo data is represented on Fig. 7 by the gray bands that show the bounds for 66 percent of 
the runs, 95 percent of the runs, and the minimum and maximum values. 

The cowl-open transient prior to the engine test resulted in the vehicle nose pitching up, whereas the nominal and 
Monte Carlo preflight analysis predicted a nose-down transient as shown in Fig. 7. The inaccurate prediction of the 
direction of the cowl-open transient was likely caused by an incorrect understanding of the unsteady aerodynamics 
caused by the cowl door opening. 

The inertial angle-of-attack response during the engine test is compared with the preflight Monte Carlo results in 
Fig. 7. During the engine test, the HXRV maintained the engine test conditions well within the requirements of 
± 0.5° angle of attack. The inertial angle of attack fell outside the Monte Carlo bounds for a few seconds during the 
engine test as the igniter was removed from the fueling profile. This event is an example of the difficulties in 
modeling the complexities of the scramjet operation, as it was not covered by the Monte Carlo preflight predictions. 
The HXRV simulation contains a simplified engine model that is essentially a table lookup based on flight condition 
and therefore, does not fully model the dynamics of scramjet operation. 

Following the scramjet experiment, the HXRV increased the angle of attack to arrest the buildup in dynamic 
pressure and heating. During the recovery maneuver, the HXRV experienced small amplitude angle of attack 
oscillations at a frequency of approximately 0.65 Hz as shown in Fig. 8. The preflight Monte Carlo predictions did 
not include any cases with sustained oscillations during the recovery maneuver. The cause of the oscillations is not 
fully understood. After the recovery maneuver, the vehicle transitioned to the descent mode. 

The vehicle maintained controlled flight and performed well during the descent. Figure 8 also shows that the 
inertial angle of attack was within the preflight predictions during the descent portion of the mission. The large 
changes in angle of attack during the descent are caused by preprogrammed PID maneuvers, which were preformed 
at integer Mach numbers. 

Figure 9 shows that the flight altitude profile was close to the nominal simulation results and was well within the 
Monte Carlo predictions. The bank angle was within Monte Carlo predictions throughout the descent except for a 
deviation late in the descent, seen in Fig. 10. At that point, the HXRV banked approximately 30° left wing down. 
The cause of this upset has not been conclusively determined, but the upset might have been caused by winds or 
unmodeled transonic aerodynamics. The impact location was well within the preflight prediction as shown in 
Fig. 11. Figure 12 shows the HXRV ground track following separation from the HXLV overlaid on the preflight 
Monte Carlo predictions. The California coastline is included in Fig. 12 for approximate scale and orientation. 

B. The Mach 10 Mission 

The HXRV generally flew within the Monte Carlo bounds for the Mach 10 flight, with several notable 
exceptions. 19 Figure 13 shows the in-flight angle-of-attack performance during the engine test. A slight pitch up, 
instead of a pitch down (as predicted by the nominal simulation), occurred when the cowl door opened. This pitch 
up occurred in both the Mach 7 and Mach 10 missions. Before the Mach 10 flight, the engine database uncertainties 
were increased to account for the discrepancy seen in the Mach 7 flight. As a result, the in-flight angle-of-attack 
performance when the cowl door opened fell within the Monte Carlo predictions for the Mach 10 flight. 

The expected longitudinal upsets, when the engine was turned on and off, were within the Monte Carlo 
predictions. During the remainder of the engine test, the HXRV maintained the required angle of attack. As seen in 
the Mach 7 flight, the inertial angle of attack fell outside of the Monte Carlo bounds during the engine test as the 
igniter was removed from the fueling profile. Figure 14 shows that during the engine test, the AMW deflections 
were near the Monte Carlo bounds, which is believed to result from an inaccurate prediction of the nominal Cmo 
(pitching moment at zero degree angle of attack) by the aerodynamic and engine databases. This discrepancy 
indicates that a larger Cmo uncertainty might have been appropriate. 

Following the engine test, the HXRV successfully closed the engine cowl door and pulled to the desired angle of 
attack as part of the recovery maneuver. After approximately 10 s, the HXRV began to slowly rock approximately 
±10° about the commanded bank angle, and the rocking grew to as much as ± 20° about the command when 
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guidance started maneuvering the vehicle in the direction of the aim point. During the same time period, the internal 
engine pressures indicated that air had started to flow through the engine, resulting in a pressure increase on the aft 
underbody of the vehicle, which pushed the nose down. The HXRV responded by commanding the elevators to 
move to a more negative deflection to bring the nose back up. Figure 15, which shows the left AMW deflection 
from flight compared to the Monte Carlo data, illustrates the primary effect of air flowing through the engine. The 
AMW deflection was outside of the Monte Carlo bounds at high Mach numbers. Another event of note during this 
time period was a step change in the elevator command of several degrees to counter a longitudinal disturbance. 
This disturbance is believed to have been caused by a venting in the fuel system and the spontaneous ignition of 
igniter in the engine during this process. The internal engine pressures dropped shortly before the Mach 8 PID 
maneuver, and the HXRV flew generally as predicted from Mach 8 and below. The AMW deflections continued to 
be outside of the Monte Carlo bounds until approximately Mach 5 when they started to track the nominal preflight 
predictions. Unexpected airflow through the engine and the unexpected spontaneous fuel ignition during the venting 
process are examples of phenomena that were not anticipated and thus not modeled as part of the Monte Carlo 
analysis. It is not known if the AMW deflections at the higher Mach numbers were outside the Monte Carlo bounds 
solely because of air flowing through the engine or if a mismatch in the predicted aerodynamics is also a cause. The 
HXRV flew as predicted during the rest of the descent after about Mach 5. 

IX. Observations and Lessons Learned 

It is fairly obvious that a vehicle simulation is only as good as the models of which it is composed. It is also 
obvious that a Monte Carlo analysis is only as good as the uncertainties used. In the case of the HXRV, some sets of 
uncertainties were known better than others. For example, the mass properties of the vehicle were measured directly 
and therefore were very well known with the uncertainty being largely driven by measurement accuracy. On the 
other hand, the HXRV was the first airframe- integrated scramjet-powered vehicle operated in flight. As a result, 
there was greater uncertainty in the aerodynamics, cowl motion transients, and engine operation because of the 
limited flight experience with this class of vehicle and the limited amount and applicability of wind tunnel testing. 
These models carried correspondingly larger uncertainties and, based on flight results, possibly should have carried 
larger uncertainties. 

The aerodynamic and engine performance uncertainties used in the HXRV simulation were primarily applied as 
a scalar to the vehicle forces and moments. The flight data indicated inaccurate predictions in the vehicle Cmo in 
several instances; therefore an improvement to the HXRV aerodynamic uncertainty model would have been the 
addition of larger biases to the total vehicle forces and moments. 

Monte Carlo analysis is very useful for predicting vehicle performance in the presence of known uncertainties. 
If, however, the uncertainties are believed to be conservative, it can provide a false sense of security and the 
perception that the vehicle will be robust to anything that could happen. Monte Carlo analysis also cannot account 
for unknown unknowns. The HXRV Mach 7 and Mach 10 missions experienced several events that were not 
anticipated and hence not modeled in the simulation. The unexpected pitch-up when the cowl door opened and the 
oscillations that occurred during the recovery maneuver during the Mach 7 mission are excellent examples of 
unanticipated events not modeled by the simulation and therefore not a part of the Monte Carlo analysis. The HXRV 
GNC team attempted to account for these unknown unknowns by maintaining conservative uncertainty bounds, 
primarily using uniform uncertainty distributions to produce a larger number of cases at the "corners," and 
performing stress tests past the uncertainty bounds. Although the exact uncertainties and unknown unknowns were 
not captured within the Monte Carlo bounds, the system was robust enough to account for them without catastrophic 
consequences. 

The Monte Carlo analysis tools turned out to be fairly complicated. Keeping track of the hundreds of model 
uncertainties and the interactions between the Matlab®-based tools and the simulation turned out to be a 
straightforward, but involved process. The simulation models and uncertainties evolved greatly between the initial 
Monte Carlo capability developed for Flight 1 and that used for Flight 3. It is extremely important that any Monte 
Carlo analysis tool be straightforward to use and easily modifiable, especially when the tool is projected to be in use 
for many years. 

The Monte Carlo analysis tools were not under configuration control until late in the project. This lack of 
configuration control was not a serious problem since modifications to the tool were typically done by only one or 
two people working closely together, however, the HXRV simulation was under configuration control from the 
beginning. It became desirable to place the tool under configuration control later in the program since it was used for 
an important preflight analysis. The majority of the uncertainty values were contained in the Matlab® scripts, which 
were a part of the Monte Carlo analysis tools, however, a number of the uncertainty values were contained in the 
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HXRV simulation and a random number was input from the tools to set the uncertainty in the simulation. There 
would have been several advantages in placing the uncertainties in the HXRV simulation. First, all of the 
uncertainties would have been in the same tool and been under configuration control from the beginning. Second, 
the Matlab”' scripts that generated each simulation run could have been much simpler if they had only to generate a 
set of random numbers for input to the simulation. An alternative approach would have been to incorporate all of the 
Monte Carlo run generation capabilities in the HXRV simulation. The main drawback to either of these approaches 
is that the simulation would have grown even more complex and potentially difficult to modify. The decision to 
place the Monte Carlo uncertainties and capabilities in the vehicle simulation is one that is specific to each project 
and depends on the complexity of both the simulation and the desired Monte Carlo analysis. 

X. Conclusion 

An overview of the Hyper-X research vehicle Monte Carlo analysis tools, model uncertainties, and analysis was 
presented. The methodology for selecting hardware validation test cases from the Monte Carlo data was described. 
The results of the preflight Monte Carlo analysis were compared with flight data, with an emphasis placed on 
highlighting modeling shortfalls. Modeling and Monte Carlo analysis lessons learned were presented. The preflight 
Monte Carlo analysis indicated that the X-43A control system was robust to the preflight uncertainties and provided 
the Hyper-X project an important indication that the vehicle would likely be successful in accomplishing the mission 
objectives. The X-43A accomplished all of the objectives during the two successful flights, however, the actual in- 
flight performance was not entirely captured within the Monte Carlo bounds. The X-43A flight control laws and 
guidance algorithms were robust enough as a result of the preflight Monte Carlo analysis, to ensure that the vehicle 
was able to handle the unexpected without undue consequences. 
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Tables 


Table 1. Monte Carlo parameters saved at 1, 10, and 100 Hz. 


Description 

1 Hz 

10 Hz 

100 Hz 

Angle of attack 


X 

X 

Sensed angle of attack 


X 

X 

Angle of attack command 


X 

X 

Estimated angle of attack 


X 

X 

Axial acceleration 


X 

X 

Side acceleration 


X 

X 

Vertical acceleration 


X 

X 

Inertial angle of sideslip 


X 

X 

Angle of sideslip 


X 

X 

Coefficient of drag 


X 


Coefficient of lift 


X 


Cowl door command 



X 

Cowl door brake 



X 

Left wing command 


X 


Right rudder command 


X 


Right wing command 


X 


Aileron command 


X 

X 

Elevator command 


X 

X 

Rudder command 


X 

X 

Right rudder deflection 


X 

X 

Left wing deflection 


X 

X 

Right wing deflection 


X 

X 

Estimated engine fuel to air ratio 



X 

Actual engine fuel to air ratio 



X 

FADS angle of attack bias 


X 

X 

FADS dynamic pressure 
estimate 


X 

X 

FADS valid flag 


X 

X 

Time since separation 



X 

Flight path angle 


X 

X 

Altitude 

X 

X 

X 

Rudder hinge moment 


X 


Left wing hinge moment 


X 


Right wing hinge moment 


X 


Reference altitude 

X 

X 


Normal acceleration command 


X 


Geodetic latitude 

X 


X 

Sensed latitude 


X 


Sensed longitude 


X 


Mach number 

X 

X 

X 

Sensed Mach number 


X 


Pressures, valve settings, etc. 



X 

Axial acceleration 



X 

Average axial acceleration 



X 

Axial acceleration at CG 



X 
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Table 1. Concluded. 


Description 

1 Hz 

10 Hz 

100 Hz 

Side acceleration 



X 

Average side acceleration 



X 

Side acceleration at CG 



X 

Normal acceleration command 


X 


Normal acceleration 


X 

X 

Average normal acceleration 



X 

Normal acceleration at CG 



X 

Roll rate 


X 

X 

Bank angle 


X 

X 

Bank angle command 


X 

X 

Sensed bank angle 

X 

X 

X 

PID complete flag 




PID start flag 




Cowl door position 



X 

Heading angle 


X 

X 

Pitch rate 


X 

X 

Dynamic pressure 

X 

X 

X 

Sensed dynamic pressure 


X 

X 

Yaw rate 


X 

X 

Range to go 

X 

X 

X 

Pitch attitude 


X 

X 

Time 

X 

X 

X 

North velocity 


X 


East velocity 


X 


Down velocity 


X 


Latitude 

X 

X 

X 

Longitude 

X 

X 

X 
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Table 2. Monte Carlo statistics from the beginning of the Mach 8 PID to splashdown. 



Run # 

Parameter 

Units 

Maximum 

Average 

Maximum 

Standard 

Maximum 

Average 

All 

Minimum 

Average 

Minimum 

Standard 

Minimum 

Maximum 

differential 

Maximum 

Minimum 

Maximum 

differential 

Angle of 
attack 

deg 

15.01 

10.3 

2.58 

6.54 

-1.48 

2.25 

1.83 

0.44 

751 

925 

596 

Angle-of- 

attack 

command 

deg 

6 

6 

0 

6 

6 

6 

0 

0 

1 

1 

1 

Estimated 
angle of 
attack 

deg 

16.67 

10.42 

2.5 

6.58 

-2.35 

2.3 

1.77 

0.44 

951 

951 

596 

Axial 

acceleration 

G 

-0.07 

-0.21 

0.05 

-0.36 

-1.22 

-0.77 

0.24 

0.03 

114 

515 

999 

Angle of 
sideslip 

deg 

3.63 

1.11 

0.42 

0.02 

-2.87 

-1.14 

0.42 

0.22 

581 

535 

951 

Left rudder 
command 

deg 

12.22 

8.89 

1.08 

0.21 

-11.96 

-9.21 

0.88 

4.17 

668 

620 

547 

Right 

rudder 

command 

deg 

12.22 

8.89 

1.08 

0.21 

-11.96 

-9.21 

0.88 

4.17 

668 

620 

547 

Left wing 
command 

deg 

20 

16.64 

3.76 

5.72 

-5.33 

-0.59 

1.54 

8.05 

2 

239 

629 

Right wing 
command 

deg 

20 

16.77 

3.64 

5.84 

-5.11 

-0.33 

1.4 

7.76 

2 

455 

629 

Aileron 

command 

deg 

7.74 

4.48 

0.94 

0.12 

-7.5 

-4.14 

1.18 

0.72 

284 

401 

785 

Elevator 

command 

deg 

23.9 

16.4 

4.45 

5.78 

-3.93 

-0.18 

1.31 

9.32 

251 

904 

629 

Rudder 

command 

deg 

8.78 

5.98 

0.91 

0.16 

-10.06 

-4.85 

0.87 

1.24 

668 

629 

951 

Right 

rudder 

position 

deg 

11.96 

8.54 

1.08 

0.06 

-11.98 

-9.12 

0.89 

1.79 

668 

620 

907 

Left wing 
position 

deg 

21.3 

16.84 

3.88 

5.99 

-4.53 

0.24 

1.6 

2.15 

546 

239 

629 
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Table 2. 


Continued. 



Run # 

Parameter 

Units 

Maximum 

Average 

Maximum 

Standard 

Maximum 

Average 

All 

Minimum 

Average 

Minimum 

Standard 

Minimum 

Maximum 

differential 

Maximum 

Minimum 

Maximum 

differential 

Right wing 
position 

deg 

21.24 

16.98 

3.77 

6.01 

-4.81 

0.27 

1.46 

2.2 

28 

904 

629 

Aileron 

position 

deg 

7.18 

4.25 

0.88 

0.03 

-7.54 

-4.15 

1.14 

0.72 

481 

907 

785 

Elevator 

position 

deg 

21.19 

16.48 

3.95 

6 

-4.26 

0.41 

1.53 

2.18 

611 

142 

629 

Rudder 

position 

deg 

11.96 

8.54 

1.08 

0.06 

-11.98 

-9.12 

0.89 

1.79 

668 

620 

907 

Right 
rudder rate 

deg/s 

1.79 

1.59 

0.06 

0 

-1.69 

-1.54 

0.06 

1.19 

907 

179 

881 

Left wing 
rate 

deg/s 

2.15 

0.62 

0.18 

0 

-1.79 

-0.64 

0.17 

1.05 

629 

999 

629 

Right wing 
rate 

deg/s 

2.2 

0.63 

0.19 

0 

-1.86 

-0.68 

0.19 

1.05 

629 

999 

25 

Flight path 
angle 

deg 

8.19 

0.15 

0.81 

-5.4 

-59.75 

-32.38 

7.58 

0.03 

90 

46 

571 

Normal 

acceleration 

command 

G 

2.5 

2 

0.18 

0.99 

-2.5 

0.36 

0.6 

5 

13 

193 

193 

Filtered 

normal 

acceleration 

G 

5.56 

1.76 

0.39 

0.96 

-0.99 

0.28 

0.33 

0.24 

951 

571 

596 

Bank angle 

deg 

100.8 

36.86 

19.94 

2.83 

-93.35 

-32.32 

11.07 

3.11 

629 

999 

629 

Pitch 

attitude 

deg 

17.08 

9.09 

2.88 

1.1 

-54.67 

-27.11 

8.19 

0.43 

90 

46 

596 

Roll rate 

deg/s 

308.43 

78.66 

42.12 

-0.01 

-247.08 

-64.02 

43.93 

13.92 

629 

999 

999 

Pitch rate 

deg/s 

43.15 

10.32 

6.19 

-0.07 

-28.14 

-10.31 

3.98 

5.38 

596 

304 

596 

Yaw rate 

deg/s 

17.96 

8.46 

2.87 

0.01 

-17.39 

-6.19 

2.91 

1.36 

386 

841 

632 

Dynamic 

pressure 

psf 

2054.93 

1371.04 

228.74 

930.64 

293.89 

677.37 

268.73 

0.46 

181 

290 

46 
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Table 2. 


Concluded. 



Run # 

Parameter 

Units 

Maximum 

Average 

Maximum 

Standard 

Maximum 

Average 

All 

Minimum 

Average 

Minimum 

Standard 

Minimum 

Maximum 

differential 

Maximum 

Minimum 

Maximum 

differential 

Sensed 

dynamic 

pressure 

psf 

1713.47 

1349.81 

202.9 

924.76 

239.26 

670.32 

270.53 

45.39 

43 

801 

43 

Estimated 

dynamic 

pressure 

psf 

3000 

3000 

0 

1136.54 

305.79 

681.92 

257.77 

11.79 

1 

526 

31 

Mach 

number 

n/a 

8.27 

7.97 

0.1 

4.38 

0.46 

0.89 

0.14 

0 

300 

90 

778 

Sensed 

Mach 

number 

n/a 

8 

8 

0 

4.38 

0.46 

0.9 

0.14 

0.01 

458 

90 

825 

Altitude 

ft 

103218 

100772 

920 

71008 

-4 

0 

1 

2 

824 

308 

299 

Angle-of- 
attack error 

deg 

3.25 

0.71 

0.55 

-0.03 

-4.76 

-1.52 

0.8 

0.09 

854 

57 

36 

Bank angle 
error 

deg 

91.86 

31.74 

10.88 

-2.79 

-99.2 

-36.23 

19.6 

4.9 

999 

629 

999 

Pitch 

attitude 

error 

deg 

53.65 

26.71 

8.04 

-1.01 

-16.73 

-8.86 

2.84 

0.54 

46 

90 

304 

Mach 

number 

error 

n/a 

0.28 

0.1 

0.05 

0 

-0.19 

-0.07 

0.06 

0 

300 

996 

951 

Dynamic 

pressure 

error 

psf 

567.83 

103.97 

103.67 

5.88 

-538.52 

-86.9 

107.39 

0.99 

181 

912 

951 
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Table 3. Monte Carlo uncertainty models and summary of parameters varied. 


Uncertainty model 

Number of parameters 

Parameters varied 

Uncertainty 

distribution 

Separation condition 

12 

Angle of attack, angle of 
sideslip, altitude, 
longitude, latitude, Mach 
number, attitude angles, 
angular rates 

Uniform 

Mass properties 

8 

Center of gravity 
position, moments of 
inertia, weight 

Uniform 

Control surface actuator 
performance 

66 

Actuator performance 
parameters for the 2 all- 
moving wings and the 
2 rudders 

Uniform 

HXRV aerodynamics 

17 

Aerodynamic 
coefficients and biases 

Gaussian 

Separation aerodynamics 

9 

Scale factors on 
separation aerodynamics 
deltas 

Uniform 

Piston performance 
(separation) 

7 

Piston performance terms 

Uniform 

Flush airdata system 
(FADS) 

8 

Pressure port CP and 
measurement uncertainty 

Uniform 

Atmospheric conditions 

1 

Month 

Uniform 

Winds 

12 

Altitude (10-120) kft 

Uniform 

12 

Horizontal speed 

Uniform 

12 

Wind direction 

Uniform 

12 

Vertical wind velocity 

Uniform 

Sensors 

51 

Gyro error and bias terms 

Uniform 

42 

Accelerometer error and 
bias terms 

Uniform 

12 

INS initialization and 
misalignment errors 

Uniform 

Engine performance 

1 

Engine performance 
scalar 

Uniform 

Propulsion performance 

4 

Propulsion performance 
parameters 

Uniform 

Total number of 
parameters 

286 
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Table 4. Monte Carlo model uncertainties. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 


Roll rate 

deg/s 

0.00 

2.00 

6.00 

-6.00 

Uniform 

Project 

requirement 

document 



Pitch rate 

deg/s 

0.00 

1.00 

3.00 

-3.00 

Uniform 

Project 

requirement 

document 



Yaw rate 

deg/s 

0.00 

1.00 

3.00 

-3.00 

Uniform 

Project 

requirement 

document 



Angle of 
attack 

deg 

0.00 

0.50 

1.50 

-1.50 

Uniform 

Project 

requirement 

document 


Separation 

condition 

Angle of 
sideslip 

deg 

0.00 

0.33 

1.00 

-1.00 

Uniform 

Project 

requirement 

document 


Pitch attitude 

deg 

UA 

UA 

UA 

UA 

Uniform 

HXLV Monte 
Carlo 



Roll attitude 

deg 

0.00 

5.00 

15.00 

-15.00 

Uniform 

Project 

requirement 

document 



Heading 

deg 

UA 

UA 

UA 

UA 

Uniform 

HXLV Monte 
Carlo 



Altitude 

ft 

UA 

UA 

UA 

UA 

Uniform 

HXLV Monte 
Carlo 



Longitude 

deg 

UA 

UA 

UA 

UA 

Uniform 

HXLV Monte 
Carlo 



Geodetic 

latitude 

deg 

UA 

UA 

UA 

UA 

Uniform 

HXLV Monte 
Carlo 



Mach 
number 
(Flight 2) 

n/a 

7.08 

0.14 

7.50 

6.65 

Uniform 

Project 

requirement 

document 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Mass 

properties 

X-Axis CG 

ft 

~6 

0.0030 

6.0089 

5.9911 

Uniform 

Measurement 
and analysis 


Y-Axis CG 

ft 

-0 

0.0022 

0.0022 

-0.0022 

Uniform 

Measurement 
and analysis 


Z-Axis CG 

ft 

-0 

0.0083 

0.0083 

-0.0083 

Uniform 

Measurement 
and analysis 


XX principle 
moment of 
inertia 

slug ft 2 

~50 

-1 

-65 

-45 

Uniform 

Measurement 
and analysis 


YY principle 
moment of 
inertia 

slug ft 2 

-800 

-15 

-850 

-750 

Uniform 

Measurement 
and analysis 


ZZ principle 
moment of 
inertia 

slug ft 2 

-800 

-15 

-900 

-750 

Uniform 

Measurement 
and analysis 


XZ moment 
of inertia 

slug ft 2 

-25 

-2 

-30 

-20 

Uniform 

Measurement 
and analysis 


X-43A 
weight at 
separation 

lb 

-3000 

-3 

-2910 

-3010 

Uniform 

Measurement 
and analysis 


Control 

surface 

actuators 

Left wing 
actuator 
natural 
frequency 

rad/s 

59.00 

8.00 

75.00 

35.00 

Uniform 

Test data and 
analysis 


Left wing 
actuator 
damping 
ratio 

n/a 

0.70 

0.14 

1.00 

0.30 

Uniform 

Test data and 
analysis 


Left wing 
actuator rate 
limit 

deg/s 

119.00 

4.00 

119.00 

107.00 

Uniform 

Test data and 
analysis 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 


Left wing 
rate limit at 
max load 

deg/s 

95.00 

2.00 

101.00 

89.00 

Uniform 

Test data and 
analysis 



Left wing 
actuator 
command 
quantization 

deg 

0.04 

0.01 

0.05 

0.02 

Uniform 

Test data and 
analysis 



Left wing 
actuator 
command 
error slope 

deg/deg 

0.00 

0.01 

0.02 

-0.02 

Uniform 

Test data and 
analysis 


Control 

surface 

actuators 

Left wing 
actuator 
command 
error offset 

deg 

0.00 

0.27 

0.80 

-0.80 

Uniform 

Test data and 
analysis 


Right wing 
actuator 
natural 
frequency 

rad/s 

59.00 

8.00 

75.00 

35.00 

Uniform 

Test data and 
analysis 



Right wing 
actuator 
damping 
ratio 

n/a 

0.70 

0.14 

1.00 

0.30 

Uniform 

Test data and 
analysis 



Right wing 
actuator rate 
limit 

deg/s 

119.00 

4.00 

119.00 

107.00 

Uniform 

Test data and 
analysis 



Right wing 
rate limit at 
max load 

deg/s 

95.00 

2.00 

101.00 

89.00 

Uniform 

Test data and 
analysis 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Control 

surface 

actuators 

Right wing 
actuator 
command 
quantization 

deg 

0.04 

0.01 

0.05 

0.02 

Uniform 

Test data and 
analysis 


Right wing 
actuator 
command 
error slope 

deg/deg 

0.00 

0.01 

0.02 

-0.02 

Uniform 

Test data and 
analysis 


Right wing 
actuator 
command 
error offset 

deg 

0.00 

0.27 

0.80 

-0.80 

Uniform 

Test data and 
analysis 


Rudder 

actuator 

natural 

frequency 

rad/s 

59.00 

8.00 

75.00 

35.00 

Uniform 

Test data and 
analysis 


Rudder 

actuator 

damping 

ratio 

n/a 

0.70 

0.14 

1.00 

0.30 

Uniform 

Test data and 
analysis 


Rudder 
actuator rate 
limit 

deg/s 

119.00 

4.00 

119.00 

107.00 

Uniform 

Test data and 
analysis 


Rudder rate 
limit at max 
load 

deg/s 

95.00 

2.00 

101.00 

89.00 

Uniform 

Test data and 
analysis 


Rudder 

actuator 

command 

quantization 

deg 

0.04 

0.01 

0.05 

0.02 

Uniform 

Test data and 
analysis 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Control surface 
actuators 

Rudder 
actuator 
command 
error slope 

deg/deg 

0.00 

0.01 

0.02 

-0.02 

Uniform 

Test data and 
analysis 


Rudder 
actuator 
command 
error offset 

deg 

0.00 

0.27 

0.80 

-0.80 

Uniform 

Test data and 
analysis 


Aerodynamics 

Axial force 
coefficient 
(CA) 

n/a 

Nominal 

(8-11)% 

+(24-33)% 

-(24-33)% 

Uniform 

Wind tunnel, 
CFD, analysis 

Varied as a 
function of 
Mach no. 

Side force 
coefficient 
(CY) 

n/a 

Nominal 

(7-10)% 

+(21-30)% 

-(21-30)% 

Uniform 

Wind tunnel, 
CFD, analysis 

Varied as a 
function of 
Mach no. 

Normal 

force 

coefficient 

(CZ) 

n/a 

Nominal 

(6-9)% 

+(18-27)% 

-(18-27)% 

Uniform 

Wind tunnel, 
CFD, analysis 

Varied as a 
function of 
Mach no. 

Roll 

moment 

coefficient 

(Cl) 

n/a 

Nominal 

(7-10)% 

+(21-30)% 

-(21-30)% 

Uniform 

Wind tunnel, 
CFD, analysis 

Varied as a 
function of 
Mach no. 

Pitch 

moment 

coefficient 

(Cm) 

n/a 

Nominal 

(6-9)% 

+(18-27)% 

-(18-27)% 

Uniform 

Wind tunnel, 
CFD, analysis 

Varied as a 
function of 
Mach no. 

Yaw 

moment 

coefficient 

(Cn) 

n/a 

Nominal 

(7-10)% 

+(21-30)% 

-(21-30)% 

Uniform 

Wind tunnel, 
CFD, analysis 

Varied as a 
function of 
Mach no. 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Aerodynamics 

Axial force 
bias 

(CA bias) 

n/a 

Nominal 

UA 

UA 

UA 

Uniform 

Wind tunnel, 
CFD, analysis 


Side force 
bias 

(CY bias) 

n/a 

Nominal 

UA 

UA 

UA 

Uniform 

Wind tunnel, 
CFD, analysis 


Normal 
force bias 
(CZ bias) 

n/a 

Nominal 

UA 

UA 

UA 

Uniform 

Wind tunnel, 
CFD, analysis 


Roll 

moment 

coefficient 

bias 

(Cl bias) 

n/a 

Nominal 

UA 

UA 

UA 

Uniform 

Wind tunnel, 
CFD, analysis 


Pitch 

moment 

coefficient 

bias 

(Cm bias) 

n/a 

Nominal 

UA 

UA 

UA 

Uniform 

Wind tunnel, 
CFD, analysis 


Yaw 

moment 

coefficient 

bias 

(Cn bias) 

n/a 

Nominal 

UA 

UA 

UA 

Uniform 

Wind tunnel, 
CFD, analysis 


Pitch 

moment due 
to pitch rate 
(Cmq) 

n/a 

Nominal 

26.67% 

+80% 

-80% 

Uniform 

Wind tunnel, 
CFD, analysis 


Yaw 

moment due 
to yaw rate 
(Cnr) 

n/a 

Nominal 

25% 

+75% 

-75% 

Uniform 

Wind tunnel, 
CFD, analysis 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Aerodynamics 

Roll moment 
due to roll 
rate (Clp) 

n/a 

Nominal 

25% 

+75% 

-75% 

Uniform 

Wind tunnel, 
CFD, analysis 


Yaw moment 
due to roll 
rate (Cnp) 

n/a 

Nominal 

66.67% 

+200% 

-200% 

Uniform 

Wind tunnel, 
CFD, analysis 


Roll moment 
due to yaw 
rate (Clr) 

n/a 

Nominal 

66.67% 

+200% 

-200% 

Uniform 

Wind tunnel, 
CFD, analysis 


Separation 

aerodynamics 

Uncertainty 
on separation 
aerodynamics 

UA 

UA 

UA 

UA 

UA 

UA 

UA 


Separation 

piston 

performance 

Uncertainty 
on separation 
piston 

performance 

UA 

UA 

UA 

UA 

UA 

UA 

UA 


Flush airdata 
sensing 
(FADS) 
system 

PPT2 

measurement 

uncertainty 

n/a 

0.00 

0.33 

1.00 

-1.00 

Uniform 

Wind tunnel, 
CFD, analysis 


PPT 3 

measurement 

uncertainty 

n/a 

0.00 

0.33 

1.00 

-1.00 

Uniform 

Wind tunnel, 
CFD, analysis 


PPT 4 

measurement 

uncertainty 

n/a 

0.00 

0.33 

1.00 

-1.00 

Uniform 

Wind tunnel, 
CFD, analysis 


PPT 5 

measurement 

uncertainty 

n/a 

0.00 

0.33 

1.00 

-1.00 

Uniform 

Wind tunnel, 
CFD, analysis 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Flush airdata 
sensing 
(FADS) 
system 

PORT 2 CP 
uncertainty 
multiplier 

n/a 

0.00 

0.33 

1.00 

-1.00 

Uniform 

Wind tunnel, 
CFD, analysis 


PORT 3 CP 
uncertainty 
multiplier 

n/a 

0.00 

0.33 

1.00 

-1.00 

Uniform 

Wind tunnel, 
CFD, analysis 


PORT 4 CP 
uncertainty 
multiplier 

n/a 

0.00 

0.33 

1.00 

-1.00 

Uniform 

Wind tunnel, 
CFD, analysis 


PORT 5 CP 
uncertainty 
multiplier 

n/a 

0.00 

0.33 

1.00 

-1.00 

Uniform 

Wind tunnel, 
CFD, analysis 


Atmosphere 

Atmosphere 
model month 

Month 

Nominal 

n/a 

12.00 

1.00 

Uniform 

Dryden 

Meteorological 

Group 


Winds 

10 kft North- 
South wind 

ft/s 

-9.44 

29.94 

80.38 

-99.26 

Uniform 

Dryden 

Meteorological 

Group 


20 kft North- 
South wind 

ft/s 

-9.92 

44.52 

123.64 

-143.47 

Uniform 

Dryden 

Meteorological 

Group 


30 kft North- 
South wind 

ft/s 

-13.16 

59.64 

165.75 

-192.07 

Uniform 

Dryden 

Meteorological 

Group 


40 kft North- 
South wind 

ft/s 

-7.68 

49.51 

140.84 

-156.20 

Uniform 

Dryden 

Meteorological 

Group 


50 kft North- 
South wind 

ft/s 

-1.21 

30.19 

89.37 

-91.79 

Uniform 

Dryden 

Meteorological 

Group 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Winds 

60 kft North- 
South wind 

ft/s 

-0.12 

16.04 

48.00 

-48.24 

Uniform 

Dryden 

Meteorological 

Group 


70 kft North- 
South wind 

ft/s 

-2.14 

9.66 

26.85 

-31.13 

Uni form 

Dryden 

Meteorological 

Group 


80 kft North- 
South wind 

ft/s 

-3.24 

9.23 

24.46 

-30.94 

Uniform 

Dryden 

Meteorological 

Group 


90 kft North- 
South wind 

ft/s 

-2.86 

11.54 

31.74 

-37.47 

Uniform 

Dryden 

Meteorological 

Group 


100 kft 
North-South 
wind 

ft/s 

-0.29 

15.94 

47.54 

-48.12 

Uniform 

Dryden 

Meteorological 

Group 


110 kft 
North-South 
wind 

ft/s 

8.07 

20.20 

68.67 

-52.53 

Uniform 

Dryden 

Meteorological 

Group 


120 kft 
North-South 
wind 

ft/s 

9.80 

24.71 

83.94 

-64.34 

Uniform 

Dryden 

Meteorological 

Group 


10 kft East- 
West wind 

ft/s 

20.53 

22.04 

86.65 

-45.60 

Uniform 

Dryden 

Meteorological 

Group 


20 kft East- 
West wind 

ft/s 

44.88 

38.63 

160.77 

-71.01 

Uniform 

Dryden 

Meteorological 

Group 


30 kft East- 
West wind 

ft/s 

67.21 

49.93 

216.99 

-82.58 

Uniform 

Dryden 

Meteorological 

Group 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Winds 

40 kft East- 
West wind 

ft/s 

83.95 

44.19 

216.53 

-48.63 

Uniform 

Dryden 

Meteorological 

Group 


50 kft East- 
West wind 

ft/s 

65.42 

26.28 

144.26 

-13.42 

Uni form 

Dryden 

Meteorological 

Group 


60 kft East- 
West wind 

ft/s 

31.54 

17.86 

85.11 

-22.03 

Uniform 

Dryden 

Meteorological 

Group 


70 kft East- 
West wind 

ft/s 

6.08 

17.85 

59.64 

-47.47 

Uniform 

Dryden 

Meteorological 

Group 


80 kft East- 
West wind 

ft/s 

-2.96 

24.45 

70.40 

-76.32 

Uniform 

Dryden 

Meteorological 

Group 


90 kft East- 
West wind 

ft/s 

3.95 

36.24 

112.66 

-104.77 

Uniform 

Dryden 

Meteorological 

Group 


100 kft East- 
West wind 

ft/s 

25.40 

49.15 

172.86 

-122.06 

Uniform 

Dryden 

Meteorological 

Group 


1 10 kft East- 
West wind 

ft/s 

60.67 

47.86 

204.24 

-82.90 

Uniform 

Dryden 

Meteorological 

Group 


120 kft East- 
West wind 

ft/s 

90.18 

56.14 

258.60 

-78.24 

Uniform 

Dryden 

Meteorological 

Group 


10 kft vertical 
wind velocity 

ft/s 

-0.02 

3.33 

9.98 

-10.02 

Uniform 

Dryden 

Meteorological 

Group 
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Table 4. Continued. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Winds 

20 kft vertical 
wind velocity 

ft/s 

-0.02 

3.33 

9.98 

-10.02 

Uniform 

Dryden 

Meteorological 

Group 


30 kft vertical 
wind velocity 

ft/s 

-0.03 

3.33 

9.97 

-10.03 

Uni form 

Dryden 

Meteorological 

Group 


40 kft vertical 
wind velocity 

ft/s 

-0.01 

3.33 

9.99 

-10.01 

Uniform 

Dryden 

Meteorological 

Group 


50 kft vertical 
wind velocity 

ft/s 

0.00 

3.33 

10.00 

-10.00 

Uniform 

Dryden 

Meteorological 

Group 


60 kft vertical 
wind velocity 

ft/s 

0.00 

3.33 

10.00 

-10.00 

Uniform 

Dryden 

Meteorological 

Group 


70 kft vertical 
wind velocity 

ft/s 

0.00 

3.33 

10.00 

-10.00 

Uniform 

Dryden 

Meteorological 

Group 


80 kft vertical 
wind velocity 

ft/s 

0.00 

3.33 

10.00 

-10.00 

Uniform 

Dryden 

Meteorological 

Group 


90 kft vertical 
wind velocity 

ft/s 

0.00 

3.33 

10.00 

-10.00 

Uniform 

Dryden 

Meteorological 

Group 


100 kft 
vertical wind 
velocity 

ft/s 

0.00 

3.33 

10.00 

-10.00 

Uniform 

Dryden 

Meteorological 

Group 


110 kft 
vertical wind 
velocity 

ft/s 

0.00 

3.33 

10.00 

-10.00 

Uniform 

Dryden 

Meteorological 

Group 
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Table 4. Concluded. 


~ denotes approximate value n/a = not applicable UA = unavailable 


Modeled 

uncertainty 

Parameter 

Units 

Mean 

1 sigma 

Maximum 

Minimum 

Uncertainty 

type 

Source 

Notes 

Winds 

120 kft 
vertical wind 
velocity 

ft/s 

0.00 

3.33 

10.00 

-10.00 

Uniform 

Dryden 

Meteorological 

Group 


Sensors 

Gyro and 
accelerometer 
performance 
uncertainties 

UA 

Proprietary 

UA 

UA 

UA 

Gaussian 

Vendor 


Engine 

performance 

Engine 

performance 

parameter 

n/a 

0.50 

0.17 

1.00 

0.00 

Uniform 

Wind tunnel, 
CFD, analysis 


Fuel system 
performance 

Fuel flow 
parameters 

UA 

UA 

UA 

UA 

UA 

UA 

Test data and 
analysis 
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Table 5. Validation test selection criteria. 

[X denotes a parameter was selected for a maximum or minimum value] 


Parameter 

Separation 

Engine Test 

Recovery 

Maneuver 

Descent 

Splashpoint 

Angle of attack 

X 

X 

X 

X 


Angle of sideslip 

X 

X 


X 


Left wing 
deflection 

X 

X 


X 


Right wing 
deflection 

X 

X 


X 


Right rudder 
deflection 

X 

X 


X 


Flight path angle 



X 

X 


Altitude 


X 




Altitude error 



X 

X 

X 

Mach 


X 


X 

X 

Normal 

acceleration 

X 

X 

X 

X 


Roll rate 

X 

X 




Pitch rate 

X 

X 




Dynamic 

pressure 


X 

X 

X 

X 

Yaw rate 

X 

X 




Phi 

X 

X 


X 


Theta 

X 


X 

X 


Time 





X 
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Figures 



Figure 1. The X-43A 3-view. 



Figure 2. Hyper-X flight trajectory. 
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MCAT 



Figure 3. Monte Carlo analysis tool overview. 
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Figure 4. Statistics for maximum filtered normal acceleration. 


Distance 

traveled, 

mi 



Drag dispersion, a 


Figure 5. Hyper-X Research Vehicle down range versus drag dispersion. 
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Figure 6. Validation test cases compared to full Monte Carlo results. 
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Figure 7. In-flight inertial angle of attack during scramjet test conditions compared to Monte 
Carlo predictions for the Mach 7 mission. 
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Figure 8. Inertial angle of attack compared to Monte Carlo predictions for the Mach 7 mission. 
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Figure 9. Inertial altitude compared to Monte Carlo predictions for the Mach 7 mission. 
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Figure 10. Bank angle compared to Monte Carlo predictions for the Mach 7 mission. 
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Figure 11. Mach 7 mission splashpoint compared to preflight Monte Carlo predictions. 
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Figure 12. Mach 7 mission ground track compared to preflight Monte Carlo predictions. 
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Figure 13. In-flight inertial angle of attack during scramjet test conditions compared to Monte 
Carlo predictions for the Mach 10 mission. 
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Figure 14. In-flight elevator position compared to Monte Carlo predictions for the Mach 10 mission. 
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Figure 15. Hyper-X Research Vehicle Mach 10 mission elevator position compared to Monte 
Carlo data. 
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